TRANSITION FROM A SUPERSONIC TO A SUBSONIC SOLAR WIND B. Durney 


ABSTRACT The transition from a supersonic to a subsonic corona was investigated by increasing the 
density at the base of the corona (initially = 9.3X10^/cm^) whfle keeping the 
temperature there constant “ 2,1X10^ °K), For the initial values of and 
the solution of the inviscid solar wind equations is of the Parker type. As iV^ is increased, 
Parker type supersonic solutions cease to exist fox NJ^ ^ 1.17X10®/cm^. No sub- 
sonic solutions exist, however, for where >3X10®/cm^. For NJ^ 

the solar wind equations permit supersonic solutions which are not of the Parker 
type, with the temperature varying as (1/r)^^^ for large distances. 


The transition region separating supersonic expansions 
of the solar wind [Parker, 1958] from subsonic 
expansions [Chamberlain, 1961] has been considered by 
Parker [1965] . Here we study this transition region for 
a given value of the temperature at the base of the 
corona. Chamberlain's notation [1961] will be used 
throughout this paper; that is, the dimensionless values 
of the temperature, square of the velocity, and the 
inverse of the radial distance will be defined as 

T = TIT I /kTi X = GMq m/kTi r ( 1) 

where Ti is a reference temperature taken equal to 
2X 10^ w the radial expansion velocity, m the average 
mass (for a hydrogen-helium mixture 10:1), Mq the mass 
of the sun, and G and k are the gravitational and 
Boltzmann constants, respectively. The usual momen- 
tum and energy equation can then be written as 
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For A we take [Chamberlain, 1961] : = 5.8X10^/C 

where C is the mass flow given in terms of N, i//, and X 
by C” X"^ is the total particle density). In 

equation (2b) epo is the residual energy per particle at 
infinity . For €oo 0 a well-known solution of equations 
(2a) and (2b) is the Parker supersonic solution. Then for 
small values of X 
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The asymptotic expansion for \]j is obtained from 
equation (2b). The parameter Dq in equation (3) is left 
undetermined; this allows equations (2) to be solved as 
follows. Given Cand epo, the integration of (2) and (3) is 

started from small values of X 
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toward the sun; the condition that the numerator and 
denominator of equation (2a) 

<t>i = \-2Tl\-dTld\ 02 = 0.5(1 -t/0) (4) 
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should vanish for the same value of X(X^) determines 
then Dp. This solution is supersonic for large distances 
and subsonic for X > \p. Therefore, given the values of C 
and €oo the integration of equations (2) determines Np 
and Tq at the base of the corona; that is, the curves 
€oo = constant and C = constant are known in the Npy Tp 
plane [Kopp, 1968] . In the present paper, C and eoo 
were adjusted to give the desired value of Tp 
(Tp = 2AX\(f°K) and increasing values, of Np 
(Np >93X\Q'^ Icm^ ). In figure 1, Dp is plotted versus 



Np\ the dashed curve represents interpolated values. In 
figure 2, Coo is plotted versus Np. Thus, from figures 1 
and 2 it is clear that Dp vanishes for 
Np =A^g = 1.17X10®/cm^, that is, for a finite value of 
Coo. It is easily seen that the vanishing of implies that 
the conductive flux at infinity is zero. Therefore, 
supersonic solutions of the Parker type cease to exist for 
>i\^. However, since ioi Np =N^ the total energy 
flux at infinity, €oo» is not zero, it is to be expected from 
physical grounds that as Np is increased the flow will 
remain supersonie. This was confirmed by numerical 
calculations. No subsonic solutions were found for 
A^^<10^/cm^. 

The above considerations led to the search for a 
supersonic solution of equations (2) with an asymptotic 
expansion different than that given by equation (3). It 
was indeed found that the full equations (2) also admit a 
solution with the following asymptotic expansions 

T = Do\^'^ [1 + X?1 X + + £>3X® + Z>4X^ 

+ AX’'^ + ...J .. . 

(5a) 

^ + 2X + Co\^'^ [1 + C,X + CjX^'® + CsX®'^ 

+ C4X* +CsX’'^ + ...] (5b) 


Figure 1. Density at the base of the corona versus Dp. 
The dotted lines are interpolated values. 
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As in the Parker solutions the parameter Dp in equation 
(5) is left undetermined, and equations (2) can again be 
solved by integrating from small values of X toward the 

0 ^ — iteration procedure was used to determine so 

^ that 0i(X) = 02 (X) = 0 (cfi eq. (4) for the same value of 

® X = X^). This determines the critical point. The iteration 

Figure 2. Residual energy at infinity versus Np for procedure is based on the fact that if 0i vanishes first 
Tp ==2.1X10^ (0i(X^) = 0 and 02(X) #O for X < X^), then Dp is too 
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X 


Figure 3. TversusX for C = 5.47X1 0'' 0.887. 



X 

Figure 4. versus X for C = 5.47X10^ 0.887. 

large and it should be decreased; if on the other hand 02 
vanishes first, then Dq is too small and it should be 
increased. This method allows the critical point X^ to be 
known as accurately as desired. Once is determined a 
Taylor expansion of r and 0 in the neighborhood of 
can be used to integrate across the critical point. In 
figures 3 and 4, r and 0 are plotted versus X for 
C = 5.47X10^ , Coo = 0.887. .With the help of the asymp- 
totic expansions (5) the integrations were started from 


X = 4.7X1 toward the sun. The arrow indicates the 
position of the sonic point. For the above values of C 
and €oo the values of the density and temperature at the 
base of the corona are Tp-l.lXlO^'^K and 
= 2.88X10^ /cm^, respectively. Therefore, to sum- 
marize, if the temperature at the base of the corona is 
assumed constant and equal to 2.1X10^'’K, the solution 
of the solar wind equations is of the Parker type for 
values of l,17XiO®/cm^ . For values of Np such 
that l.llXlO^jcm^ <N<NSiNS>3Xm the 

flow continues to be supersonic but differs from the 
Parker solution. The asymptotic expansions of the 
velocity and temperature are now given by equations 
(5). For N>N^ the flow becomes subsonic. The value 
of was not determined with great accuracy; the 
numerical calculations suggest that N^'^IXIO^ jcm^ . 
Since the integrations are time consuming no supersonic 
solutions were evaluated for Np > 3X 10® /cm® . 

As the density at the base of the corona is increased 
from an initial value of = 9.3X 10^/cm® and the 
temperature there is held constant (7^ = 2.1X 10^°K) 
there is not a direct transition from a Parker type 
supersonic flow to a subsonic flow. As Np is increased 
the conductive flux at infinity e£j decreases sharply and 
vanishes for a finite value of the total energy per particle 
at infinity eoo. This is not surprising, since disNp increases 
the mass flow C also increases, and since Tp has been 
kept constant Coo must decrease. The sharp decrease in 
arises as more and more of the conductive flux at 
infinity is converted into kinetic energy of the particles. 
The Parker type supersonic solutions cease to exist when 
vanishes. The total energy flux at infinity remaining, 
however, different from zero, it is to be expected from 
physical grounds that if Np increases still further there 
will not be an immediate transition to a subsonic flow 
but rather a transition to a different type of supersonic 
flow (cf. eqs. (5) and (6)). The expansion becomes 
adiabatic [Chamberlain, 1961] and the heat conduction 
flux is much smaller than the internal energy flux 
[Hundhausen, 1971], Any further increase in Np must 
now be entirely accounted for by a decrease in the 
kinetic energy flux at infinity. The velocity reaches a 
maximum and decreases towards a finite value as r ^ 
Finally, for values of Np large enough, Coo vanishes and 
the flow becomes subsonic. 
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